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We consider a Bianchi type / model in which anisotropic dark 
O I energy is interacting with dark matter and anisotropic radiation. With 

this scenario, we investigate the vahdity of the generahzed second 
law of thermodynamics. It is concluded that the validity of this law 
I depends on different parameters like shear, skewness and equation of 

^ ■ state. 
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, , 1 Introduction 

X 

' The observational evidence that our universe is making a transition from a 

decelerating phase to an accelerating is the major development in cosmology. 
Supernova la data [H |2] gave the first indication of the accelerated expansion 
of the universe. This was confirmed by the observations of anisotropies in 
the cosmic microwave background (CMB) radiations as seen in the data from 
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satellites such as Wilkinson Microwave Anisotropy Probe (WMAP) [3] and 
large scale structure [4] . Astrophysical observations indicate that accelerated 
expansion of the universe is driven by an exotic energy with large negative 
pressure which is known as dark energy (DE). 

The data indicates that the universe is spatially flat and is dominated 
by 76% DE and 24% other matter (20% dark matter and 4% other cosmic 
matter). Despite all lines of observational evidence, the nature of DE is still 
a challenging problem in theoretical physics. Several models have been pro- 
posed such as, quintessence [5], phantom field [6], tachyon field j^, quintom 
[8] and the interacting DE models, Chaplygin gas [9], holographic models 
[TU] and braneworld models [TT], etc. However, none of these models can be 
regarded as being entirely convincing so far. 

In black hole physics, the temperature and entropy are proportional to 
the surface gravity at the horizon and the area of the horizon [121 US] re- 
spectively. Also, the temperature, entropy and mass of the black hole satisfy 
the first law of thermodynamics [H]. This leads to the relationship between 
the black hole thermodynamics and the Einstein field equations. Jacobson 
[15j derived the field equations from the first law of thermodynamics for all 
local Rindler causal horizons. Padmanavan [16] formulated the first law of 
thermodynamics on the horizon using the field equations for a general static 
spherically symmetric spacetime. Later, this relationship was developed in 
the cosmological context by taking universe as a thermodynamical system 
bounded by the apparent horizon. Cai and Kim [TTj showed the equiva- 
lence of the first law of thermodynamics to the Friedmann equations with 
the generalized second law of thermodynamics (GSLT) at the horizon. 

Cosmological models in which DE interacts with dark matter and other 
cosmic matter are well-known in the literature [TH] . The discovery of black 
hole thermodynamics has led to the thermodynamics of cosmological models. 
Bekenstein |l3j proved that there is a relation between an event horizon and 
thermodynamics of black hole. The event horizon of black hole is in fact the 
measure of entropy which is generalized to the cosmological models so that 
each horizon corresponds to an entropy. The GSLT is generalized in such a 
way that the sum of the time derivative of each entropy must be increasing. 
People studied the thermodynamics of different models in which different 
cosmic constituents interact. 

Wang et al. |19j found that both the first and the second law of thermo- 
dynamics break down at the event horizon. The first law of thermodynamics 
is restricted to nearby states of local thermodynamic equilibrium while the 
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event horizon reflects global features of the spacetime. Also, due to the exis- 
tence of the cosmological event horizon, the universe should be non-static in 
nature and, as a result, the usual definition of the thermodynamical quan- 
tities on the event horizon may not be as simple as in the static spacetime. 
Finally, they have argued that as the event horizon is larger than the appar- 
ent horizon so the universe bounded by the event horizon is not a Bekenstein 
system. 

There is a large body of literature [20]- [32] dealing with the thermody- 
namics of the universe bounded by the apparent horizon and the validity of 
GSLT. Mazumder et al. [33], [31] have explored the validity of GSLT using 
the first law of thermodynamics and found some conditions on its validity. 
Debnath [53] investigated the validity of GSLT in the scenario of Friedmann 
Robertson Walker (FRW) model in which DE interacts with dark matter 
and radiation. Mubasher et al. [36j proved the validity of GSLT for all time, 
independent of geometry and equation of state (EoS) parameter. In a recent 
paper [37], the validity of GSLT is explored in the Kaluza-Klein cosmology 
with modified holographic DE. 

Jacobs [38] studied the spatially homogeneous and anisotropic Bianchi 
type / (BI) cosmological model with expansion and shear but without ro- 
tation. He discussed anisotropy in the temperature of CMB and cosmic ex- 
pansion both with and without magnetic field. Akarsu and Kilinc [39] inves- 
tigated anisotropic BI models in the presence of perfect fluid and minimally 
interacting DE with anisotropic EoS parameter. They found that anisotropy 
of the DE did not always promote the anisotropy in the expansion. 

In this paper, we take the BI model in which anisotropic DE interacts 
with isotropic dark matter as well as anisotropic radiation and investigate 
the validity of GSLT. The plan of this paper is the following: in next section, 
we take the BI model as the representation of the universe and evaluate 
shear parameters {R,S), skew parameters {6, 7, (, ^). Also, we formulate 
the corresponding field equations, the generalized Friedmann equation and 
equations of continuity in the effective field theory. Section 3 presents the 
analysis of the validity of GSLT. The last section concludes the paper results. 
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2 Anisotropic DE Interacting with Two Flu- 
ids 

The metric representation of the BI universe is given as 

ds^ = -de + A^{t)dx^ + B^{t)dy^ + C\t)dz\ (2.1) 

where A, B, C are scale factors along the x, y, z axes respectively. When 
A = B = C , the BI model reduces to the flat FRW model. Thus BI is 
the generalization of the fiat FRW model. The energy-momentum tensor is 
defined as 

= Tl' , + T/;, + r/', , (2.2) 

V {m)u iX)'^ ^ ' 

where the subscripts m, A, x denote dark matter, anisotropic DE and 
anisotropic radiation, respectively, and corresponding energy momentum ten- 
sors are given as follows: 

'P{m)v ~ ( PraiPmiPm^Pm)-, (2-3) 
P[k)v = {-pA,PWx,P{A)y,P{A)z), (2.4) 

'^{x> = i-Px^Pix}^^P{x)y^Pix}z)- (2.5) 
The corresponding EoSs for these components are 

(2.6) 
(2.7) 
(2.8) 

To parameterize the anisotropy of pressure, we define skewness parame- 
ters. These parameters actually measure the deviation in pressure along the 
y and z-axis from the x-axis. The skewness parameters for anisotropic DE 
are defined as [10] 

^ P{A)x - P{A)y 2^ ^ P(A)z - P{A)x 

Pa ' Pa ' 

Using Eqs.( l2.6p and ( 12. 9p . we can write Eq.f l2.4p as follows: 

P{A)iy = [-PA:^APA, (c^A + 35)pA, (^A + 37)pa]. (2.10) 



P{A)x 


— PA^(A)x-, 


PWv - 


PA^{A)y, 


PiA> = 


'- PA^{A)z, 


P(m)x 


Pm^{m)x) 


P(in)y = 


~ Pm^{m)y 1 


P(m)z 


Pm^(m)z) 


P{x)x 


= Px^ix)^^ 


P{x)y = 


Px^{x)y^ 


P{x)z = 


Px^(x>- 
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The skewness parameters for anisotropic radiation are 

Px ' Px ' 

With the help of Eqs.dM]) and f l2lT|) . Eq.([22D turns out to be 

^(x)^ = [-Px'^xPx> K + 3C)Px. K + 30Px]- (2-12) 
Making use of Eqs.(Ein]) and fl^TT^ in Eq.(|23), we have 

Tl^ = [- Pm, t^mPm, t^mPm, t^mPm] + ["PA, ^APA, (^A + 36)Pa, {ua + ^1)Pa\ 

+ [-Px.^xPx.K + 3C)Px>K + 30Px]- (2-13) 
The corresponding field equations take the following form: 

Ab BC CA ^ ^ , 

m^BC^CA = ^P- + P- + P^^^ (2.14) 

B C BC , , , , 

B^C^M^ = -[w^Pm + WAPA + WxPx]' (2-15) 

A c Ac 

A^C^T^ = -[w,np,n + (a;A + 35)pA + (wx + 3C)Px]' (2-16) 

A B Ab 

A^B^AB ^ -knPm + (c^A + 37)pA + (^^x + 30Px]- (2-17) 
The continuity equations, T^'^-i, = 0, are 

PA + Hpa[3{1 + ua) + 5(3 -2R + S)+ 7(8 -2S + R)] = -Q', (2.18) 

p^ + 3Hp^{l + ojm) = Q. (2.19) 

Here Q and Q' are introduced because two principal components of the uni- 
verse are mutually interacting, which lead to some loss in other cosmic compo- 
nent. We take this component as radiation for which the continuity equation 
becomes 

p^ + Hp^[3{l + oj^) + C(3 -2R + S)+ e(3 - 2S + R)] = Q' -Q. (2.20) 
Here, we define the mean expansion rate as an average Hubble rate H by 



The difference of expansion rates as the Hubble normahzed shear parameters 
R and S are given as 



I ,A B, ^ I A C , , , 

Shear parameters are also deviation measuring parameters which parameter- 
ize anisotropy in scale factors by taking expansion along standard. 
Here, all expansion rates are positive and R > —3, S" < 3. This is also true 
either for R = or S = 0. The generalized Friedmann equation takes the 
form [ID] 

tt2 _ Pm + PA+ Px /r, 

- 3(1_ l(i?2 + 52_^^))- ^^-^^^ 

We take [12] Q' = FaPa, Q = T„iPk- The ratios for energy densities as 

Pm Py 

ri = — , ra = — . 
Pa Pa 

The corresponding EoS parameters in the effective field theory will become 
_ ^ , Ta 5i3-2R + S) + ji3-2S + R) 

= ^m-^, (2.25) 
3Hri 



Consequently, the continuity equations fl2.18l) - fl2.20p turn out to be 



PA + 3HpA{l + ul") = 0, (2.27) 
p^ + 3HpUl + u;'£) = 0, (2.28) 
p^ + 3ifp^(l + <^) = 0. (2.29) 



3 Generalized Second Law Of Thermodynam- 
ics 

Now we investigate the validity of the GSLT in BI universe bounded by 
apparent horizon with size L which coincides with Hubble horizon in the 
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case of flat geometry, i.e., L = j^. The first law of thermodynamics gives 

PdV + dE 



TdS = PdV + dE, dS 



T 



where T, S, E and P are the temperature, entropy, internal energy and 
pressure of the system respectively. The corresponding entropies will become 

_ P^dV + dEA _ P„4V + dEm ,„ _ P^dV + dE^ 

"•JA — j: ) dbm — J, , ~ ' 

(3.30) 

where Fa, Px^ Pm, -Ea, E^ and Em are the pressures and internal energies of 
anisotropic DE, anisotropic radiation and dark matter, respectively. We as- 
sume that the system is in equilibrium, which implies that all the components 
of the system have the same temperature given by [13] 



T = ^, (3.31) 

where K is the surface gravity of the black hole. At equilibrium, the horizon 
and other components of the system have same temperature. Thermody- 
namical quantities are related to the cosmological quantities by the following 
relations: 

Pa = oo'/^pA, Px = ujfpx, Pm = Uj'^^Pm, (3.32) 
Ea = ABCpA, E^ = ABCp^, Em = ABCpm, (3.33) 

where V = ABC is the volume of the system containing all the matter. 

The entropy of the horizon is Sh = where A* is the surface area of 
the black hole and k is the Boltzmann's constant. Here A* = AttL'^ so that 
we have Sh = krcL"^, which leads to 

Sh = 2k7rLL. (3.34) 

Also, the time derivative of Eq. fl3.30p yields 

• _PaVj^ . P^V + 4 . PmV + Em 

Using Eqs.dMlD-lESnD, ( 1334]) and ( l335|) . it follows that 

Stotai = 2knLL, (3.36) 
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where Stotai is the sum of the matter entropy and horizon entropy. Now 
making use of Eqs.(E2ZD-(I22SD and ^TM . we obtain 

• H _L^Ha L^H^[R{2R-S) + S{2S-R)] 

18^ ' ^^-^^^ 

where 

a = (1 + Ua)pa + (1 + Wm)Pm + (1 + ^x)Px 

+ 1[{PaS + PxC)(3 -2R + S) + (pA7 + PxO(3 - 2S + R)], (3.38) 

and 

R^ + S'^ — RS , . 

13 = 1 > 0. (3.39) 

Inserting this value of L in Eq.( l3.36p . it follows that 

. kTcHL^a k'KH^L^[R{2R- S) + S{2S - R)] . ^ , 

Stotai - . (3.40) 

This implies that if Stotai ^ 0, then GSLT holds. Thus the validity of GSLT 
depends upon the skewness, shear and state parameters. From Eg. (13. 40 p . 
Stotai ^ implies that 

knHL^a _ knH^L^[R{2R - S) + S{2S - R)] ^ 

After substituting the values of a and /3 in the above equation, it follows that 

9(p + p)> [R{2R -S) + S{2S - R)] - 3[{paS + p^0{3 -2R + S) 
+ (PA7 + PxO(3-25 + i?)]. (3.41) 

Now, we consider the following three interesting cases involving the condition 
on different parameters: 

I. ^ = 0, II. /3 > 0, III. R = S = 6 = -r = ^ = C = 0. 



Case I 

Here H and Stotai tend to infinity for /? = from Eqs. fl2.23p and (13.401) 
respectively. It may happen for very large time (t — ?■ oo), i.e., when the 
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expansion rate is very high. In this case, all the useable energy in the universe 
will be converted into another form of energy which is not useable. This 
stage is also known as the heat death of the system. The heat death of the 
universe is said to be a suggested fate of the universe. At this time, all the 
thermodynamic free energy will be diminished from the universe and motion 
or life cannot sustain any more. In the language of physics, the entropy in 
the universe will reach to its maximum value. 

Case II 

In this case, by takeing /3 > in Eq. fl3.4Up . we observe that the validity of 
this law depends on the skewness, shear and EoS parameters. Further, if 
we remove anisotropy in the expansion, its validity depends on skewness as 
well as equation of state parameters, and if we remove anisotropy in fluid, 
then it depends on shear and equation of state parameters. Hence, GSLT is 
conditionally valid. 

Case III 

When we take R = S = 5 = -f = C = C = 0,'Bl reduces to the fiat FRW 
model. Using these values in Eq. (13.411) . we have {p + p) > 0, which is the 
null energy condition in the FRW universe. Consequently, Eq. fl3.40p leads 
to Stotai = kT^ip + p) and hence Stotai > for {p + p) > 0. Thus we can 
say that GSLT holds for all time in this case, if the null energy condition 
for the considered matter is satisfied in the FRW universe. This case implies 
that the validity of GSLT depends on the null energy condition which is not 
discussed in [SB]- Therefore using results produced in this case, we can say 
that the validity of GSLT depends on null energy condition and background 
geometry. 

4 Conclusion 

This paper is devoted to study the validity of GSLT with the BI universe 
model. We assume that anisotropic DE is interacting with dark matter and 
anisotropic radiation. Using this scenario, we have formulated the conditions 
under which GSLT is valid. We have discussed three cases for its validity. 
For a particular value of /3 = 0, the first case provides the general validity 
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of GSLT. In this case, we see that at infinite expansion the heat death of 
universe will take place. At this stage all types of motion and life will be 
finished. The second case /? > gives conditional validity of this law. In 
the third case, BI reduces to the fiat FRW metric and the general validity 
of GSLT depends upon the condition {p + p) > 0, which is the null energy 
condition. We would like to mention here that the validity of this law for the 
FRW universe with the same scenario has been proved [SB] independent of 
background geometry, EoS parameter and interacting scenario. This analysis 
indicates that the validity of GLST for FRW depends on geometry as well as 
null energy condition. We conclude that GSLT is conditionally valid in the 
BI type universe with anisotropic perfect fiuid. 
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